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r— n ABSTRACT 

Oh 

Taking advantage of an improved genetic optimization algorithm for fitting 
unconstrained Kepler orbits to the Hipparcos intermediate astrometric data, we 

1 obtain additional orbital solutions for 81 Hipparcos stars with previous stochas- 
tic (failed) solutions. The sample includes astrophysically interesting objects, 
including the old disk wide binary HIP 754, the nearby AGB star HIP 34922 (L2 

^ Pup), and the nearby M2 dwarf HIP 5496 (GJ 54, at 8 pc from the Sun) which 

has a resolved M dwarf companion. 

O Subject headings: astrometry — binaries: general 

o 

> 

^ 1. Introduction 

A genetic optimization-based algorithm of unconstrained orbital solutions on Hipparcos 



Intermediate Astrometry Data (HIAD, ESA 1997) was developed and applied to 1561 stars 



from the annex of stochastic solutions, resulting in 65 orbiting binaries at the 99% confidence 



level (Goldin & Makarov 2006). Most of the 65 objects had not been known as binaries, 
although at least two stars on the list of "new" binaries turned out to be parts of a long- 
term spectroscopic monitoring carried out by other researches (R. Griffin, J. Sperauskas, 
priv. comm.). Since our first publication about orbiting stochastic stars, we have realized 
that a simple improvement of the algorithm can be effected if the intermediate astrometric 
data produced by the two data analysis consortia (NDAC and FAST) are combined prior 
to running the optimization program, as opposed to our initial decision to run the program 
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separately for the NDAC and FAST data points and to compare the final results. Typically, 
each great circle measurement produced a pair of abscissa data points, one derived by FAST, 
and the other by NDAC. Our initial disinclination to average the data derived from the same 
observations was justified by the strong correlation which is underestimated in the catalog. 
At the same time, some of the astrometric abscissae derived from the same observations 
differ significantly between the two consortia, betraying unknown model or systematic errors. 
These errors can be diminished by simple averaging, but the estimated precision can only 
slightly improve because of the correlated noise. The safest and most conservative approach 
is to compute the mean values of the common data points, but to adopt the larger of the 
two formal errors as the expected standard deviation of the result, effectively assuming a 
correlation of 1. Thus, our choice of methods and algorithms is driven by the intent to 
minimize the number of false positives (which can not be precluded completely with these 
noisy data) and to produce confidence estimates as realistic as possible. 



Uncertainties and bias 



Our results are affected by stochastic uncertainties, as well as by considerable bias, es- 
pecially in the important parameters clq (apparent orbit size) and e (eccentricity). Therefore, 
it is of crucial importance not only to compute an orbital fit to a set of data points, but also 
to evaluate the probability distributions for the derived parameters. The one-dimensional 
nature of Hipparcos measurements, the strong irregularity of observing times and position 
angles, and the presence of unmodeled disturbances, such as visual companions, all combine 
to often produce oddly shaped, asymmetric distribution functions. In order to compute 



realistic precision and bias estimates, we use the parametric bootstrap method (Efron & 



Tibshirani 1993). 



Our general optimization method is described in (Goldin & Makarov 2006). First, we 



minimize x on abscissae residuals given observed transit times tj and transit times error 
estimates o^. The best fit yields estimated times t ei such that t\ = t J ei + a^ej, where ej are 
normally distributed normal variables with variance 1 and j is the number of realization. 
With sufficiently large number of realizations, we build a histogram for each estimated fitting 
parameter. The expected error of the estimate can be derived from the width of the his- 
togram. In case of zero bias and nearly symmetric distribution, the mode of the distribution, 
its mean and median are close to value derived from the original data. 

As an example, we consider the eccentricity fits produced by our unconstrained op- 
timization algorithm for a subset of 231 stars identified as orbital solutions in the main 
Hipparcos catalog (ESA 1997). Obviously, for high values of eccentricity the bias on the 
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positive side is limited by the hard limit of 0.99 imposed by our program. But if we consider 
only solutions with small enough eccentricity (less then 0.5), we can clearly see that only 
one solution out of 90 has a negative estimated bias. Only 3 solutions out of these 90 with 
small eccentricity have a negative estimated bias in semimajor axis ao- 

Due to the nonlinear mapping of t{ to the orbital parameters, the distributions of the 
mass function and ao may have extremely long tails. Let V est be the estimated parameter 
value from the original data, and V 3 - the estimate from a bootstrap realization j. To avoid 
the unwanted influence of occasional extreme values of V on our error and bias estimate, we 
define bias via the median, viz., B(V) = V es t — median(V), and la-error bars (which can be 
asymmetric) as cr_(V) = median(U) — qo.\e(V) and cr + (V) = qo.siiV) ~ median(U). where 
q g (V) is the g-quantile of the distribution of V J . 

For sufficiently small errors and sufficiently large number of data points, it is possible to 
estimate the unbiased value of V un uased = V est — B{V). We do not attempt this correction for 
the stochastic solutions considered in this work. For the 231 stars with catalogued Hipparcos 
orbital solutions, one can see from Fig. [T] that the bias of e is consistently positive if e est 
is not too large. Note that our eccentricity estimation is restricted to the interval [0, 0.99], 
but the bias estimation is restricted only at the upper limit. It is likely that this bias is 
responsible for the abundance of solution with eccentricities close to the cutoff value of 0.99. 



3. Orbits 

An automated run of the modified optimization algorithm on all 1561 stochastic stars 
in the Hipparcos catalog resulted in about 200 new solutions with estimated confidence 
levels greater than 99%. The meaning of this confidence estimation is that the HIAD data 
for these stars is perturbed in such a way that the apparent positions are consistent with 
the specified orbital model, with the probability of a random occurrence of such data for 
a single, unperturbed star less than 1%. This is not to say that the given solutions are 
accurate to 99%, or even that the stars are 99%-probability binaries. Indeed, we find that 
roughly one half of all automated solutions from the stochastic set have eccentricities equal 
to 0.99, and large orbital sizes. This unexpected set of failed fits indicates that undetected 
orbital motion in our somewhat idealized model only accounts for a fraction of stochastic 



solutions. Fabricius & Makarov (2000) demonstrated that a number of stochastic stars 
are visual, resolved binaries with negligible orbital motion, where the secondary companion 
was not listed in the Hipparcos input catalog, or given at an incorrect position. Using 
correct initial assumptions about the relative position of components in binary and multiple 
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Fig. 1. — Orbital eccentricity and bias of eccentricity estimated by the parametric bootstrap 
method for 231 binary stars with orbital solutions in the Hipparcos catalog. 
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systems can yield regular, high-quality astrometric solutions by a straightforward least- 
squares adjustment. But the model of fixed visual companions is another simplification, 
bound to meet limited success with the real objects presenting a host of complicating features. 
For example, a visual companion can have a detectable motion around the primary, resulting 
in an effectively different proper motion; the inner pair in a hierarchical triple system can 
be orbiting; one or both components of an orbiting pair can be variable stars evoking the 



variability-imposed motion (VIM) effect on the observed photocenter (Pourbaix 2001). For 
red stars, the data itself may need recalibration, since the assumed color was in many cases 
biased ( Platais et al.|2000~ ). In principle, increasingly sophisticated algorithms can be created 



to deal with such cases; in practice though, the optimization problem becomes cumbersome 
with the available data because we are quickly running out of data points with each additional 
nonlinear fitting parameter. 

At this point, we simply rejected all the spurious-looking solutions with extreme eccen- 
tricities. The remaining 81 fits are specified in Table [TJ Orbital fits for the stars published in 
our preceding paper dGoldin fc Makarov |[2006| are also omitted, even if they differ for some 



stars. The estimated astrometric parameters, period P in days, eccentricity e, periastron 
time T in days, longitude of the ascending node u in degrees, position angle of the node 
Q in degrees, and inclination i (zero for face-on orbits) in degrees are specified in the table. 
These data are based on 1000 genetic optimization trials for each object. Confidence limits 
for each fitting parameter are computed at several confidence levels from the set of 1000 
simulations, but only the 84th and 16th percentiles are used to compute the robust uncer- 
tainty limits in Table [TJ which correspond to the ±ler interval of N(m,a). These intervals 
provide a minimalistic estimation of the parameter uncertainties, but carry no information 
about the shape of the distribution, which can be strongly non-Gaussian, as discussed in the 
previous paragraph. Our solutions are perfectly unconstrained, and are not based on any 
prior information except the HIAD. 

A literature search for the stars listed in Table 1 reveals that three objects have been 
known as spectroscopic binaries from precision radial velocity observations. Two of them 



have accurate spectroscopic orbits published, viz., HIP 34164 = G 108-53 (Latham et al. 



2002] and HIP 67480 = e Boo ( |Griffin|[l985| . Both stars were also re- analyzed by Jancart 



et al. (2005), whose computations were also based on the HIAD, but the essential nonlinear 
parameters (P, e and uj) were fixed at their spectroscopic values. The cited spectroscopic 
data, our unconstrained fits and constrained astrometric solutions from ( Jancart et al.||2005 ) 



are specified in Table [2] for comparison. Encouragingly, our periods come out fairly close 
(within 1-2 a) to the more accurate spectroscopic estimates, and the eccentricities are even 
closer to the mark for both stars. The unconstrained estimates of uj probably include the 
uncertainty of 180° in both cases. Therefore, we do not see any reasons to be concerned with 
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our method and results from the comparison with the spectroscopic data. The constrained 



solutions for the apparent orbit size a and inclination i from (Jancart et al. 2005) are 
much stronger than ours, because fixing the troublesome nonlinear parameters (especially, e) 
simplifies the multidimensional x 2 surface and removes many of the valley-like local minima. 
For HIP 34164, our estimate of % is comfortably close to the constrained parameter. The 
corrected mass function, Mf/(Mi + M 2 ) 2 , is fairly high at 0.095 M for this high proper 
motion, moderately metal-poor ([Fe/H]= —0.68) system. 

For the star HIP 67480 (e Boo), our unconstrained fit and the constrained solution from 
( Jancart et al.||2005 ) produce similar orbit sizes, but significantly different inclination angles. 



The former value of % is brought into agreement with the small radial velocity amplitude 



reported by Griffin (1985) for this K4III star and the spectroscopic aisinz w 14.1 Gm. 
Our unconstrained solution would imply a companion of moderate mass in the absence of a 
spectroscopic data, whereas we are probably dealing with an orbit seen almost face-on. 



Finally, the F7V star HIP 110785 = HD 212754 has been listed by |Griffin fc Suchkov 



(2003) as a single- lined spectroscopic binary, but without orbital solution. This star has 
been drawing considerable attention of observers since it was identified as a likely member 
of the Hyades moving group. Our solution yields an almost edge-on orbit of size a$ « 0.43 
AU. From this, and the period P = 878 d, we conclude that the mass ration M2/M1 may be 
fairly small (~ 1/5). 



The Geneva-Copenhagen (GC) spectroscopic survey of Hipparcos solar-type stars (Nord- 



strom et al. 2004) provides an independent verification of our results. Most of the stars in 
this survey obtained only a few radial velocity (RV) measurements, which was insufficient to 
derive an orbital fit, but sufficient to establish the fact of spectroscopic binarity. For stars 
with more than one measurement the probability is specified that the RV has not changed 
significantly relative to the estimated measurement error. These probability estimates are 
useful to make sure that most of our orbital solutions indeed correspond to binary systems, 
although the quality of our solutions can not be verified. We find 15 stars in common be- 
tween Table 1 and GC with estimated probabilities of constant RV, viz., HIP 754, 1768, 
3645, 12726, 17482, 32307, 49638, 56447, 64790, 84696, 101430, 110785, 113177, and 113699. 
All these probabilities are zero or very small in GC. We conclude that the incidence of 
spectroscopic binarity among the new solutions is probably high. 



The stars HIP 3645, 43067, 49638, 50567, 84696, 109095 are listed in (Makarov & Kaplan 



2005 ) as astrometric binaries whose long-term proper motions from the Tycho-2 catalog are 
statistically different from the short-term Hipparcos proper motion. Such discrepancies are 
usually caused by significant long-period orbital motion. The fitted periods for these stars 
are indeed at the upper end of the range in Table 1, exceeding 3-4 years. Admittedly, this 
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is not a clean check on the quality and reliability of our solutions, because the stochastic 
proper motions in Hipparcos can be perturbed more than their formal errors manifest. 

The data in Table 1 indicate that our method is especially sensitive for orbital periods 
between 1 and 4 years, with relatively few shorter or longer periods detected. The smallest 
orbital size appears to be 5 mas. It is difficult to derive robust orbital solutions for smaller 
apparent orbits even for the brightest stars, where the typical single measurement precision is 
comparable to or less than 5 mas. This is an interesting fact pertaining to the scarcity of new 
low-mass companions detected from the HIAD, and the complete lack of planet signatures. 
Rather than a hint of lower than expected precision of the data, we consider it to be a 
still poorly understood feature of the given nonlinear optimization problem, which has been 
empirically found to have a worse condition than what could be naively expected from the 



number of data points and the number of fitting parameters (see also Pourbaix 2001). 



The orbital fits for a few stars in Table 1 are so definite and robust, that one may wonder 
why they were not resolved by the Hipparcos data analysis team and placed in the annex of 
orbital solutions. Fig. [2] shows the reconstructed orbit and the averaged observations for one 
of these stars, HIP 113699. Our fit yields an orbital period of nearly 2 years and an apparent 
motion of a$ = 9.5 mas for this poorly investigated GO dwarf. The star was part of the 
systematic attempt to resolve Hipparcos stochastic binaries with speckle interferometry by 



Mason et al. (2001), which was unsuccessful. The estimated orbit size in combination with 
the parallax (23 mas) indicates that the companion may be much fainter than the primary. 
The estimated inclination suggests that the star should be easily detectable as SB1 in RV 



measurements. 



4. Objects of note 



HIP 5496 = GJ 54 is a nearby M dwarf at a distance of ~ 8 pc from the Sun. 
The stochastic Hipparcos parallax 123 ± 8 mas is somewhat uncertain, but a recent ground- 



based observation puts the star confidently within 10 pc (Henry et al. 2006). Golimowski 



et al. (2004) resolved a visual companion to this star using HST NICMOS observations 
and found a magnitude difference of ~ 1 mag between the components. Our solution for 
this star is relatively robust, with a period of 427(+9 — 8) d and apparent semimajor axis 
ao = 19.0(+7.4, —2.1) mas. The small orbital size for this nearby star seems to suggest a low- 
mass companion, perhaps a brown dwarf, raising doubts that the orbiting companion and the 



visual are the same one. However, the provisional data from (Henry et al. 2006) confirm our 
orbital solution, estimating the period at ~ 1.1 year and the apparent orbit of the photocenter 
at 16 mas. Obviously, the companion is bright enough (magnitude difference of 1.04 in V) to 
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Fig. 2. — The apparent orbit of the photocenter of the new binary HIP 113699. The ob- 
served abscissa residuals are shown with straight segments connecting the estimated orbital 
positions (indicated with black circles) and the actual observations. The mean epoch posi- 
tions of the Hipparcos grid of slits is marked with short straight segments. The x 2 -adjusted 
parameters (slightly different from those specified in Table 1) are: P = 667 d, e = 0.1, 
T = 216 d, % = 36.4°, a = 8.7 mas, u = 97°. 
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pull the photocenter closer to the center of mass, making the observed ellipse much smaller 
than the true one. A simple calculation shows that for a plausible total mass of 0.5 M , the 
full orbit is a = 0.88 AU, and since the observed photocenter orbit is 0.15 AU, the primary's 
orbit is 0.40 AU, and the mass ratio q = MijM\ = 0.83. Thus, the Hipparcos measurements 
and our orbital solution are consistent with the known properties of the resolved companion. 



What remains puzzling is the overall brightness of the system. According to (|Henry et 

9.82, I 



al.| |2006| ), the magnitudes are V = 9.82, I = 7.32, J = 6.00, and K = 5.13 mag. This 
compares well with the accurate photometry from (Koen et al. 2002), V = 9.812 ± 0.008 



and I = 7.357 ± 0.009 mag. Fig. [3] depicts the well-defined empirical main sequence of 
nearby dwarfs in the Mk vs. V — K axes from (Henry et al. 2004) (lower curve). The short 



curled curves with beads on them indicate the joint photometry of unresolved main-sequence 
binaries on a grid of primary absolute magnitudes, for magnitude differences between the 
companions of AV^ = +oo, 6, 3, 2, 1, and mag. The latter value corresponds to the case 
of two equally bright companions (twins). With the trigonometric parallax determined by 
Hipparcos, the system is much brighter in K than the upper envelope of possible magnitudes 



for binary stars. The parallax measured by Henry et al. (2006) (141.2 ± 3.4 mas) brings it 
closer to the range of normal magnitudes, but interestingly, our updated astrometric solution 



for parallax is close to the original Hipparcos value. Henry et al. (2006) contend that this 



system represents almost the worst case, when the orbital period is close to one year, and it 
becomes difficult to distinguish between the parallactic and orbital ellipses. The covariance 



analysis technique on linearized condition equations described in ( Goldin & Makarov 2006 ) 



is useful to estimate the degree of possible perturbation of the parallax solution. We compute 
a significant correlation coefficient of —0.24 for the pair of parameters ao and n. It can be 
suspected that the trigonometric parallax is underestimated; therefore, the orbit size may 
be slightly overestimated. We also find, surprisingly, large positive correlations between n 
and the angles u and Q of 0.57 and 0.61, respectively. It follows, that the actual difficulty 
with this star is that the parallax determination is strongly coupled with the orientation of 
the orbit, which is poorly constrained by our data. 

HIP 34922 = L 2 Pup is probably the nearest asymptotic giant branch (AGB) star to 
the Sun. AGB stars are in the terminal stages of stellar evolution characterized by massive 
loss of mass. By virtue of its large-scale, stable variability of brightness, the star has been 
classified as semiregular pulsating variable. Our unconstrained fit for this star produced a 
robust solution with a period of 141(+2, —1) days, which is identical to the variability period 
of 140.6 d. This prompted us to investigate if the star could be a Variability-Imposed-Motion 
(VIM) binary, whose photocenters move along the line connecting the components in phase 
with the variability, because of the shift of the center of light. A general VIM fit produced a 
meagre reduction in the residual x 2 ; ruling out a simple VIM model. Another possibility of 
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Fig. 3. — The empirical main sequence of field dwarfs from ( Henry et al.||2004| ) (lower curve) 
and the location of binaries with joint photometric magnitudes shown with short curves and 
dots for AV = (upper end points), 1, 2, 3, 6 and +oo mag. The loci of the star GJ 54 are 
indicated with open circles for two parallaxes, 123 and 141 mas. 
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a spurious astrometric fit is predicated by the incorrect V — I color indeces that had been 
assumed for some of the red stars in the data reduction. This could generate significant 



errors in color-dependent calibration terms, as explained in detail in (Platais et al. 2000). 
For this star, however, the assumed color (V — I = 4.18) was not too different from the 
color derived from subsequent photometric observations (V — I = 3.51), and the calibration 
error is not expected to exceed 1 mas, which is much smaller than the fitted orbit (9.5 ± 



[+4.1 — 1.3]). Thus, L 2 Pup is quite likely a real orbiting binary. Jura et al. (2002) obtained 
and studied mid-infrared images of L 2 Pup and found an extended asymmetric morphology 
of the shell. The asymmetry was later confirmed from the interferometric observations by 



Ireland et al. (2004). Lebzelter et al. (2005) detected periodic variation of radial velocity, 



but unfortunately, the small number of data points does not permit to estimate the phase 
shift between the light curve and the RV curve. We propose that the stable period of the 
variability curve, the recent dimming of the object and other peculiar properties of the 
object can be explained by a bright companion orbiting the primary AGB star and getting 
periodically obscured by its thick dust shell. A more detailed discussion of this system will 
be published elsewhere. 

HIP 754 = HD 471 has been netted in earlier searches for field subdwarfs because of its 



considerable proper motion (e.g., Sandage & Kowal 1986). However, the star is moderately 



metal poor at [Fe/H]= —0.16 (Nordstrom et al. 2004), has a small Galactic eccentricity of 



0.11, and its maximum excursion from the Galactic plane is limited to 258 pc (Allen et al. 



2000). Thus, it is not a subdwarf, but Allen et al. argue that it is one of the oldest known 



wide binaries with a photometrically estimated age of 11.6 Gyr. The wide companion is 
separated by a little less than 30". This is definitely a physical pair, with proper motions from 
the extensive Nomad catalog (Zacharias et al. 2003) (fj, a *,^s) = (174.9, —146.7) ± (1.8, 1.3) 



mas yr" 1 for HIP 754 A and (171.0,-147.5) ± (0.6,0.9) mas yr" 1 for HIP 754 B. The 
wide companion has a separation of at least 1200 AU, and a period > 34000 yr. Dynamical 
evolution and the origin of old wide binaries are interesting issues motivating the long-term 
effort to increase the sample of such systems and to collect more observations. Given the 
considerable attention to this star in the literature, it is surprising that binarity of the 
primary in this pair has remained unnoticed. Our solution suggests a fairly large apparent 
orbit of 9.8 ± (+3.9, —1.5) mas and a nearly edge-on position, explaining the large scatter 



of radial velocities observed in (Nordstrom et al. 2004). The corrected parallax is 24.8 ± 1.0 



mas. Assuming a total mass of 1.5 M , the mass of the primary is roughly 1.1 M & and the 



secondary 0.4 M Q . Tokovinin et al. (2006) argue that new stars do not form in short-period 
binaries, and we mostly owe the presence of spectroscopic binaries with periods less than 10 
d to the dynamical interaction between the inner pair and the remote tertiary companion in 
hierarchical multiple systems, which takes giga-years to shrink the inner orbit, and perhaps, 
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to eventually merge the close pair. The HIP 754 system, if it really is very old, appears to 
be an example to the contrary, when the tertiary failed for some reasons to shrink the close 
inner binary. 

Facility: HIPPARCOS 

The research described in this paper was in part carried out at the Jet Propulsion Lab- 
oratory, California Institute of Technology, under a contract with the National Aeronautics 
and Space Administration. This research has made use of the SIMBAD database, operated 
at CDS, Strasbourg, France. 
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Table 2. Comparison of astrometric fits with known spectroscopic orbits and constrained 

solutions. 



HIP 34164 



Element 
Spectroscopic 
P [d] 

e 

fM 
u> [°] 

a\ sin i [Gm] 
Astrometric, unconstrained (this work) 
P [d] 643jlS 
0.27±° ;I? 

i [°] 104±if 
,+65 
-53 
,+19.8 

-9 ^ 



HIP 67480 



ao [mas] 



612.3 ± 3.0 
0.273 ±0.015 
0.0826 ±0.0061 
248.9 ± 3.7 
91.9 ± 2.4 



111- 1 
9.3 4 



944 ± 8 
0.41 ± 0.09 
0.00013 ± 0.00005 
359 ± 15 
14.1 ± 1.9 



+ 153 

89 
+0.35 
0.21 
13 
15 

111+ 45 
1±1 -42 
+6.0 
1.8 



1167 
0.38 
134 



9.6 1 



Astrometric, constrained (Jancart et al. {2005} 

i[°] 107.4 ±8.5 174.0 ±0.5 

a [mas] 8.77 ± 0.96 7.3 ± 0.9 



Note. — Spectroscopic data from jLatham et al.| 
[2002l > and ( |Griffin|1985] > 



